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The light curve of the semiregular variable L 2 Puppis: I. A recent 
dimming event from dust 
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ABSTRACT 

The nearby Mira-like variable L2 Pup is shown to be undergoing an unprecedented dim- 
ming episode. The stability of the period rules out intrinsic changes to the star, leaving dust 
formation along the line of sight as the most likely explanation. Episodic dust obscuration 
events are fairly common in carbon stars but have not been seen in oxygen-rich stars. We 
also present a 10-/im spectrum, taken with the Japanese IRTS satellite, showing strong sil- 
icate emission which can be fitted with a detached, thin dust shell, containing silicates and 
corundum. 
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> . 1 INTRODUCTION 



L 2 Puppis (HR 2748; HIP 34922) is a bright nearby red giant with 
a pulsation period of about 140 d. Its spectral type of M5eIII and 
luminosity of 1500 L© indicate that it is evolving towards the tip of 
the Asymptotic Giant Bra nch (AGB). Evid ence for mass loss at a 
rate of 3 x 10 -7 M yr _1 (tfura et al. 2002b supports this. L 2 Pup is 



possibly the nearest star in this evolutionary phase, at a Hipparcos 
distance of 61 ± 5pc. Among known long-period AGB stars, only 
R Doradus has a similar distance. At 12 microns, L2 Pup is among 
the 15 brightest sources in the IRAS point source catalogue. 

L2 Pup is unusual in several respects. Firstly, it shows a high 
degree of optical polarization, with a variable wavelength depen- 
dence that implies a long timescale for th e growth and dissipatio n 
of dust grains (of the order of a decade; Magalhaes et al. 1986 ). 



Secondly, CO measurements by Kerschbaum & Olofsson (1999) 
indicate a very low expansion velocity (about 2.5 km s - ), which 
led them to label L2 Pup as an extreme case, with one of the small- 
est expansion velocitie s ever measured for an AGB star. The slow 
wind from L2 Pup led Winters et al. (200C, 2002) to suggest that 
this star could represent their B-model, in which mass loss is driven 
entirely by pulsations, without any significant input from radiation 
pressure on^itis^j^aim^^lris^Jias^c^n^Jto 
et al. (2002), who modelled the mass loss and suggested that the 
pulsations may be non-radial. 
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Thirdly, as we report here, this star has shown a remarkable 
change in mean visual magnitude over the past century, and is cur- 
rently undergoing a dramatic dimming. We present visual and in- 
frared photometry which characterizes this behaviour, and argue 
that the most likely cause is the formation of dust along the line of 
sight. We also present the first 10-micron spectrum of L2 Pup, ob- 
tained with the Japanese IRTS satellite, which shows strong silicate 
emission. 



2 LIGHT CURVE 
2.1 Visual observations 

We have analysed visual observations of L2 Pup from the fol- 
lowing sources: the Royal Astronomical Society of New Zealand 
(RASNZ; 12100 measurements by 110 observers, including 1100 
by A. Jones), the Variable Star Observers League in Japan (VSOLJ; 
792 measurements by 6 observers) and the Association Francaise 
des Observateurs d'Etoiles Variables (AFOEV; 389 measurements 
by 3 observers). Only data from individual observers contributing 
30 or more observations were used and we did not attempt to cor- 
rect for offsets between observers. The top panel of Fig. |j| shows 
the combined data, binned to 10-day averages. 

The variation in mean magnitude is evident. Before 1945 the 
data are too patchy and/or noisy to reach strong conclusions. How- 
ever, it is clear that the star faded significantly after 1960, recover- 
ing around 1975 but remaining fainter by about 0.5 mag compared 
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Figure 1. Wavelet analysis for L2 Pup. Top to bottom: the light curve (in 10-day bins), the WWZ transform, the semi-amplitude and the period. 



to the 1950s. Recently, starting around 1994, a dramatic fading has 
occurred. Even at maximum, the star now remains fainter than mag- 
nitude 6, and it is no longer a naked-eye variable. 

We have used wavelet analysis to search for period and am- 
plitude changes. Wavelets have been used prev iously to study 
lovig-^mo^^Mi^o^^^^^^^^^^^^^^^^^&AdiTig et al. 
1998; |Kiss et al. 1999| ). We used the weighted wavelet Z-transform 



(WWZ; Foster 1996), developed specifically for unevenly sampled 



data. We experimented with different values for the parameter c, 
which de fines the trade off between time resolution and frequency 
resolution Foster (1996), and settled on c = 0.005 as a good com- 



promise. More details of the appli cation of the WWZ tr ansform to 
long-period variables are given by Bedding et al. (1998). 



The lower three panels of Figure [l] show the wavelet plots for 
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L2 Pup, based on the light curve in the top panel. The second panel 
shows the WWZ transform, with the grey scale indicating the sig- 
nificana^of^MicJilranjM et al. 

1998). Only a small range of frequencies is shown - there is no evi- 
dence for significant power outside this range. The third and fourth 
panels show, for each time bin, the semi-amplitude (in magnitudes) 
and period (in days) corresponding to the peak of the WWZ in the 
second panel. 

It is clear from Fig. ^ that the period of L2 Pup has remained 
within a fairly narrow range (135-145 d) for the past 75 years. Such 
period jitter is common among semiregular and Mira variables. 
There is a slight indication of a lengthening of period before the 
onset of the recent dimming, from 136 to 142 days, but the change 
(at most 4%) is well within the range of normal Mira variables. 
The amplitude, on the other hand, has changed significantly over 
the years, but the changes are apparently uncorrelated with the pe- 
riod jitter. 

L2 Pup is classified as an SRb variable, a class of stars with 
poorly defined periodicity. In fact, the period of this star is remark- 
ably stable. The characteristics of the wavelet plot show that L2 Pup 
should be classified as SRa, a class closely related to Miras, often 
differing only in having smaller amplitudes. 



2.2 Early observations 

L2 Pup was di scovered to be variable by Gould in 1 872 ( Cannon & 
Pickering 1907). |Roberts (1893[ ) estimated aperiod of 137.2 d from 
the dates of nine maxima observed by himself . Gould and Williams 
during 1872-92. Cannon & Pickering (1907 1, in the Second Cata- 
logue of Variable Sta rs, quoted a period of 140.1 5 d, which they 
attributed to Roberts. |Cannon & Pickering (1909 ) listed dates of 
maxima going back to 1872, and also adopted a period of 140.15 d. 
It therefore seems that the period of L2 Pup has remained within a 
narrow range since its discovery. 

The vi sual magnitude a t maximum during the late 1 800s was 



Table 1. Near-infrared photometry of L2 Pup 



about 3.6 (Williams 1 897 ; Cannon & Pickering 1909), consis 



tent with modern (1950s) out-of-decline values. M tiller & Hartwig 
(1922) listed the magnitude range since discovery as 3.4^1.6 at 
maximum and 5.8-6.2 at minimum. This suggests that, between 
1872 and 1918, L2 Pup was never as faint as it currently is. 



2.3 Near-infrared photometry 

We have^t^akeji^fn^Ljrrfrare^jto et al. 

(2000), together with more recent observations obtained with the 
same system (see Table [j]). The infrared observations are much less 
frequent than the visual ones and do not sample the pulsation cycle 
very well. Nevertheless, the mean infrared magnitudes can be es- 
timated, and we have shown these in Fig. ^ by fitting a low-order 
polynomial to each waveband. 

The dimming at V is followed at all infrared bands, although 
with reduced amplitude. The J ~K colour has reddened during the 
dimming. 

The mean K magnitude bef ore the dimming was —2.34 (Ta- 
ble 4 of Whitelock et al._200C ). More recently, the magnitude 
dropped to —2.0 or fainte^J. The Hipparcos distance allows com- 
parison with the A'-band period-lum inosity (P-L) relation for 
nearby Miras ( Whitelock & Feast 2000} : 
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Table 2. Pre-dimming and present mean magnitudes and colours of L2 Pup 





pre-dimming 


present 


V 


4.8 


6.8 


J 


-1.03 


-0.25 


H 


-1.94 


-1.25 


K 


-2.34 


-1.8 


L 


-2.91 


-2.7 


V - K 


7.0 


8.5 


J ~K 


1.31 


1.55 


™bol 


0.73 a 


1.3 a 



assuming no extinction correction 



M K = -3.47 log P + 0.84. 



(1) 



For L2 Pup, this relation predicts Mk = —6.64, which compares 
relatively poorly to the observed value of Mk = —6.16. The recent 
fading has worsened the agreement. Thus, L2 Pup appears to be one 
of a small fraction of stars located below the P-L relation. We will 



return to this in Sec. 3.3. 

Table ^ shows the mean colours and magnitudes before and af- 
ter the dim ming. The pre-dimming NIR magnitudes are taken from 
Table 4 of Whitelock et al. (2000). The current values are derived 



1 [Bedding & Zijlstra (199g| ) gave K ■ 
error and should have read —2.15. 



-2.65, but this was a transcription 



from the fits in Fig. g, and are somewhat more uncertain because of 
the poor sampling of the pulsation cycle. 



3 DISCUSSION OF THE DIMMING 

What causes the large variations in mean magnitude of L2 Pup, 
especially the current dramatic fading? We consider three possibil- 
ities: changes in temperature, luminosity and extinction. 



3.1 Change in stellar temperature and/or luminosity 

In late M-type giants, the optical flux is very sensitive to temper- 
ature. Even a slight cooling would lead to a sharp drop in V via 
the formation of temperature-sensitive molecules (TiO, VO) that 
absorb strongly at optical wavelengths (Reid & Goldston 2002). 

We first need to estimate the effective temperature of L2 Pup. 
In the absence of a measured angular diameter, this requires refer- 
ence to published temperature scales. The temperature scales for 
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Figure 2. Infrared and visual photometry for L2 Pup. 



M-giants differ greatly for stars with and without extended at- 
mospheres. Assuming no extended atmosphere, the pre-dimming 
V — K cojou^jnd^x^o^^^u^h^pji^^f^ ^ 3200 K (B essell 



et al. 1998 
also give a J 



Ho udashelt et al. 2000^ 7 Models with this temperature 
K colour consistent with the observed value. Jura 



et al. (2002) adopted a temperature of 3400 K, also taken from pho- 
tospheric model s, which is also consistent with the scale of van 
Belle et al. (1996)! 

However, stars with extended atmospheres have m uch lower 
colour temperatures ( Feast 1996 ; van Belle et al. 2002 ). The fol- 
lowing relation was given by ^east(T99*6 ) for Miras, based on in- 
terferometric angular diameters: 



log T off = -0.474(J - K) + 4.059. 



(2) 



For L2 Pup, the pre-dimming J — K colour index yields the much 
lower value of T e g = 2775 K. Since this star is Mira-like and 
shows evidence for mass loss, the atmosphere is likely to be ex- 
tended and the lower temperature is probably a better description 
of its optical/NIR spectral energy distribution. 

During the fading, J—K increased to 1.55, which corresponds 



to a temperature decline of 660 K on the scale of Feast (1996). The 
increase in V — K to 8.5, would correspond t o a temperature de- 
cline of about 300 K in the Bessell et al. (1998) models. The bolo- 
metric correction, BCx, increased from 3.05 to 3.21 on the low- 
tempe rature scale, and from 3.05 to 3.15 on the high temperature 
scale (IWhitelock et al. 20001 iBessell et al. 1998b. The luminosity 



decline in either case is about a factor of 1.8. Using L cx i? 2 T c * ff , 
we conclude that the radius would have decreased by 10% on the 
high-T scale and increased by 20% on the low-T scale. 

Such large changes in radius can be ruled out because of the 
stable period. The maximum period change allowed by the obser- 
vations is about 6% (see Fig. |l|). Since the period of a pulsating 
star scales inversely with the square root of stellar density, we can 
infer that the radius cannot have changed by more than 4% (since 

R 3/2 ). In fact, we can take the relative stability of the pulsa- 
tion period over the past 130 years to be evidence that the density of 
the star has not changed significantly. The period variations that are 
observed in L2 Pup are similar to those seen in many semiregular 
variables and are most likely due to the influence on the pulsation 
driving by random convective excitation (Christensen-Dalsgaard 
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et al. 200 1). We discuss this in more detail in Paper II (B edding 
et al. 2002). 

The period is stable on the same timescales over which the 
mean magnitude is variable, which argues that radius changes can- 
not be part of the explanation of the dimming. If we assume a con- 
stant radius, we could still produce the observed change in m^ \ 
if the temperature decreased by about 400 K. Such a temperature 
change cannot be ruled out but would require fine tuning. Further- 
more, such a large change in stellar luminosity over a time scale 
of a decade only occurs in AGB stars in the immediate aftermath 
of a thermal pulse. Finding any star in such a short-lived phase is 
very unlikely. In any case, the thermal-pulse scenario can also be 
ruled out because, during thi s phase, the radius strongly contracts 
and the period shortens (e.g. Vassiliadis & Wood 1993 ), which is 
not observed in L2 Pup. 

It is therefore unlikely that either the radius or the luminosity 
of the star has changed. This points towards variable extinction as 
the cause of the dimming. 



3.2 Extinction from circumstellar dust 

The plausible cause of the fading in L2 Pup is extinction by cir- 
cumstellar dust, which would not affect the amplitude or period of 
pulsation of the underlying star. Extinction laws for interstellar dust 
are well studied and predict less extinction as one moves to longer 
wavelengths. This is certainly the case for L2 Pup, as can be seen 
from Fig. [| the decline is much less dramatic in JHKL than in V. 

Direct evidence for extinction comes from the present J — K 
and K — L, whi ch are among the reddes t of the M-type variables 
in the sample of Whitelock et al. (200C). All stars with such red 



colours have K-bund amplitudes larger than 0.7 mag, while the 
amplitude for L2 Pup is only 0.29 mag at K. In fact, both AK 
and SHp (0.71) are small compared to M-type variables. This large 
discrepancy between colour and amplitude suggests that the red 
colours are partly caused by extinction. 

L2 Pup is also anomalously red for its period. The average 
relation between colour and period for Miras is 



J - K = -0.39 ± 0.15 + (0.71 ± 0.06) logP, 



(3) 



( Whitelock et al. 200( ). For L2 Pup, this relation predicts J — K 



1.13, compared with 1.31 pre-dimming and 1.55 at present. The 
observed values suggest that significant extinction existed even be- 
fore the present dimming. Assuming a linear relation between the 
J — K excess and the visual extinction, and assuming that the post- 
1995 dimming of 2 mag at V and 0.24 in J — K is due to extinction, 
the pre-dimming J — K excess of 0.18 would correspond to a visual 
extinction of about 1.5 mag at V. 

This estimate suggests that the long-term V-band variations 
may be due to variable extinction. The brightest epoch, in the 
1950s, could in this case give the best indication of the unreddened 
magnitude. However, while L2 Pup is Mira-like in some properties, 
it is certainly not a true Mira and so results based on applying Eq. ^ 
should be treated with some caution. 



3.3 Extinction curve and the Mira P-L relation 

Typical exti nction constants for interstellar dust are listed in Bessell 
et al. (1998). For the observed post-1995 J — K reddening, they 
give AAk = 0.15. This is less than half the decline at K shown in 
Fig. ^, suggesting this extinction curve is too steep. To investigate 
this quantitatively, we show in Fig. [3] the locus (M, V) of the mean 
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Figure 3. Infrared versus V photometry for L2 Pup, using the polynomial 
fits shown in Fig. ^ The vertical axis is the infrared magnitude, for M = 
J, H, K, L. Written on each plot are t he observed slope an d the expected 
slopes, the latter taken from Table 3 of ( ^ardelli et al. 1989[ ). 




Figure 4. Infrared colour-colour diagram of L2 Pup, using the polynomial 
fits shown in Fig. ^| The diagonal line shows the interstellar reddening vec- 
tor, with A v = 1.0, E(J -K) = 0.164 and E(K - L) = 0.045. 



infrared magnitude versus the V magnitude, for M — J, H, K, L. 
Also shown are stra ight lines with slopes taken from Table 3 of 
Cardelli et al. (1989), which are standard values for extinction by 



interstellar dust. We can see that the observed slopes are in all 
cases larger than expected, and the discrepancy is progressively 
greater at longer wavelengths. The interstellar extinction curve ap- 
pears inappropriate: instead the extinction appears to be 'greyer.' 
This could be due to a contribution from oxides, which can form 



at high temperatures (Henning 1996). The infrared colour-colour 
diagram, shown in Fig. ^, confirms that the reddening is less steep 
than expected from interstellar extinction. 

The colour differential suggests that Ak/Av ~ 0.3. If we 
assume that the pre-dimming (1995) extinction at V was Ay w 
1.5 mag, the if-band would already have been affected by 0.5 mag, 
so that Mk,o ~ —6.7. This is close to the value expected from the 
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Mira P-L relation, and well within the 1-tr uncertainty contributed 
by the error on the Hipparcos parallax. It therefore appears possi- 
ble that the Mira P-L relation is an accurate indicator of the de- 
reddened magnitude of L2 Pup. 

3.4 Dust location 

For a pulsation velocity of 3 km s~ 1 , the gas travels at most 1R, per 
year. The short timescale of the dimming thus suggests that the new 
dust is located within a few stellar radii o f the star, possibly less . 
The rapid variation in polarization angle (Magalhaes et al. 1986) 



also suggests that dust forms close to the star. The timescale for dust 
formation derived by these authors based on the polarization (about 
a decade) is consistent with the timescale of extinction variations. 

The polarization shows that the dust is not distributed isotrop- 
ically, and the variation in polarization angle shows that the distri- 
bution is not constant. Dust may form continuously, but in different 
directions: the present dimming event would, in that case, corre- 
spond to dust formation along the line of sight. 




Wavelength (/^m) 

Figure 5. The mid-infrared spectrum taken by the MIRS on board of the 
Japanese IRTS satellite. Open circles indicate channels affected by glitches. 
The error bars give the 1 -a scatter of the signal in the off source region. 
Other sources of flux uncertainty, for example, possible change of detector 
sensitivity, are not included. 



3.5 Extinction episodes in other Mira variables 

Dimming episodes lasting a number of pulsation cycles (decades) 
are seen in carbon Mi ras, with similari t ies to the dust obsc uration 
events of R CrB stars ( Percy et al. 1990 ). Lloyd Evans (1997| ) found 
that these episodes are accompanied by strong emission in the Swan 
C2 band 5 165 A, while Whitelock et al. (1997) found them to occur 



in carbon stars with thick circumstellar envelopes. 

Could th e fading of the non- carbon star L2 Pup be a related 
phenomenon? Lloyd Evans (1997) listed five Miras with dust fad- 
ings which probably have relat ively low C/O ratios (W Aql, U Cyg, 
Y Del, S Lyr and RU Vir), and |Whitelock et al. (2000| ) added Y Vel 
to this list. Unlike L2 Pup, all these stars have long periods (433- 
490 d). And, in any case, inspection of the visual light curves of 
all six stars on the AAVSO Web site do not show any evidence for 
dimmings of the magnitude and longevity of those seen in L2 Pup. 

Dimming episodes are seen in sy mbiotic Miras, wh ere they 
are attributed to variable dust extinction (Whitelock 1987). A typi- 



cal dimming at J is 1-2 mag, even larger than that seen in L2 Pup. 
We conclude that among non-symbiotic oxygen-rich Miras and 
semiregulars, the dimming of L2 Pup appears to be unique. 



4 MASS LOSS 



Whitelock et al. (2000) noted that L2 Pup stands out as one of only 
two low-amplitude SRs with a large K — [12] excess (the other 
being V CVn), implying a larger mass loss than is usual for SRs. 
There may be an effect from the low expansion velocity: the dust 
will remain close to the star longer, thus increasing the [12]-^tm 
excess. 

Stronger evidence for on-going mass loss would come from 
a 10-fim silicate feature. Such a feature has not been reported, be- 
cause of the lack of a published 10-micron spectrum. The lack of 
an original IRAS LRS spectrum is com mon for bright stars (e.g., 
o Ceti). However, IVolk & Cohen (19891) failed to recover the IRAS 



spectrum, which seems to have suffered from a technical glitch. 
ISO did not observe this star. 

L2 Pup was included in the recent survey by the Japanese 
Infra red Telescope in Space (IRTS) satellite (M urakami et al. 
1996). Obser vations were made using the Mid-Infrared Spectrom- 
eter (MIRS) ( Roelligetal. 1994 ; Yamamura et al. 1996 ), which has 



an 8 x 8 arcmin 2 aperture and scanned with a speed of 4 arcmin/sec. 
Spectra were taken on 1995 March 30 and April 5, i.e., just before 
the dimming event. We re-performed the data reduction to correct 
for effects caused by the high brightness of the star. The coverage 
in this part of the sky by the IRTS was also below average, resulting 
in only two well-centred scans over L2 Pup. The averaged spectrum 
is shown in Fig. |5| 

The spectrum shows a surprisingly strong silicate feature, 
with a band-to-continuum ratio of about 2.5. Thi s give a sil- 
icate strength parameter ( Schutte & Tielens 1989] ) of Aiq = 
2.5 log^obs/Jcont) ~ 1. The continuum between 6 and 8 /im 
has a slope consistent with the stellar temperature but is brighter 
than expected from the L-band magnitude of —2.9. 

There is no evidence for a large dust halo around L2 Pup. We 
inspected a 4 square degree field in t he four IRAS images made 
using the IRAS Software Telescope (Assendorp et al. 1995): the 



field shows complicated cirrus structure but no large envelope. 



4.1 Mass-loss rate 

The silicate feature is very strong compared to the limits on the 
underlying dust continuum. To quantify this, we fitted the spec- 
trum using the dust model of Siebenmorgen & Kruegel (1992). The 



model uses as input a blackbody, which we assigned a temperature 
of T off = 2900 K and luminosity L, = 2.0 x 10 3 L Q . The dust was 
described as a r~ 2 distribution (i.e., constant mass-loss rate) with 
an inner and outer radius: the mass-loss rate was converted to den- 
sity assuming an expansion velocity v cx _ p — 2.5 kms . The dust 
was included as silicate grains with a size distribution n oc a -3,5 , 
with diameter a > 150 A. In the fitting procedure, we adjusted the 
stellar temperature and luminosity, the dust density (or mass-loss 
rate) and inner radius, to reproduce the MIRS spectrum. The stellar 
parameters are largely determined by the shorter wavelength obser- 
vations. 

The 25 and 60-/im DIRBE measurements are much higher 
than the corresponding IRAS values. Part of this can be related to 
the different response curves: both these DIRBE pass bands cut off 
a few microns bluer than the IRAS bands. However, the increase 
over the predicted IRAS flux for the model spectrum is only about 
10%. Although the DIRBE beam is much larger than that of IRAS 
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Table 3. Fitted parameters to the infrared spectrum (1995: pre-dimming) 
of L2 Pup, corresponding to the right panel of Fig. ^. 



Fn.e 

F12 

F 25 

F 60 

F100 



Model 

(Jy) 
2120 

1810 

2100 

1260 

127 

22 



IRAS 

(Jy) 



1990 

783 

88 

27 



DIRBE 

(Jy) 



2137 
1385 
151 



Keck (flura et al. 2002|) 

(Jy) 

2500 
1700 



Parameters: 

u 



Ml 



cm 
cm 



2900 K 
2000 Lq 

7.2 x 10 14 

1.3 x 10 15 
350-200 K 
300A 

5 x lO" 10 M yr _1 
5.3 X 10- 7 M Q yr- 1 



a range between smallest grains at inner radius and largest grains at outer 
radius 

b assuming [Fe/H]= -0.7 



(as the telescope was only 19cm, compared to 60cm for IRAS), but 
there are no other other bright point sources in the beam. 

It is difficult to reconcile the photometry with the silicate fea- 
ture. This is illustrated in Fig. |^, where the left panel shows a fit 
based the silicate feature, and the right panel shows a fit that in- 
cludes all photometry. The two models differ in their outer radius, 
which is smaller for the right panel. The fitted parameters of the 
right panel (our preferred model) are listed in Table ^. The large 
number of free parameters means that the fits are not unique: i.e., a 
larger mass-loss rate with a thinner shell could also fit the spectrum. 

The most constraining aspects of the spectrum are the lack of 
dust emission at 8/^m and the strong silicate feature, which preclude 
the presence of hot dust. In our 'best' model, the inner radius is 
at a dust temperature of 350 K, well below the dust condensation 
temperature. The 17/im flux (also dominated by silicate emission) 
is also surprisingly sensitive to the inner radius and agrees with a 
dust shell detached from its condensation point. The outer radius 
is not as well determined as the inner radius. However, a single 
thin shell can fit the long-wavelength IRAS/DIRBE flux densities. 
There is no evidence for a cold distant dust shell, consistent with 



the low CO mass-loss rate found by Olofsson et al. (2002) 



Between 4 and 8/im, the spectral energy distribution indicates 
a colour tempera ture consistent with the photospheric temperature, 
lura et al. (2002) argued that dust emission dominates from 5 /im, 



based on evidence that the variability behaviour is very different 
longward of 5 /im. However, the emission in this band in Mira vari- 
ables is dominated by water vapour at ~2000 K, located in the ex- 



tended envelope at about 1.5 stellar radii (Matsuura et al. 2002) 



This water vapour also raises the flux in this region above that ex- 
trapolated from the photospheric JHK magnitudes, as seen in the 
fit. Our findings do not indicate a substantial contribution from dust 
to the spectrum shortward of 8 fjm. 

We derive a (pre-dimming) mass-loss rate of 5 x 
10~ 7 Movr~ 1 . This is similar to that estimated by J ura et al. 
(2002) (post-dimming). They found a dust mass-loss rate 4 times 
larger than derived by us because they adopted a dust outflow ve- 
locity of lOkms -1 rather than 2.5 kms -1 , but used a much lower 



gas-to-dust ratio, which almost cancels this difference. Our high 
gas-to-dust ratio is based on the assumption that L2 Pup is a thick- 
disk star (see Sec. i.3), for which we have assumed [Fe/H]= —0.7. 



In reality the gas mass-loss rate is not well determined. 



4.2 Silicate and Corundum 

The model shows that the photometric data imply a narrower sil- 
icate feature than is observed. The most likely explanation is that 
th e observed 10- ^tm feature is not purely silicate. In the scheme of 



Speck et al. (200C ), this feature in L2 Pup would be classified ei- 
ther as 'broad' or possibly as 'silicate D'. (The difference is subtle 
and requires continuum subtraction, for which we lack sufficient 
wavelength coverage.) They fitted the broad feature with a mixture 
of olivine (MgFeSiO^O and porous amorphous alumina (or corun- 
dum) (AI2O3). Corundum has the highest condensation tempera- 
ture of any circumstellar solid, and is the first condensate to form at 
T ~ 1600 K. Silicate s may form subsequently from g as-solid re- 
actions with alumina (Helens 1990). frrlron et al. (1997 ) suggested 
that the broad feature is seen in less evolved Mira variables. This 
may reflect the possiblity that, at low M, the conversion of alumina 
to silicate (replacing the Al with Mg) remains incomplete, due to a 
low reaction rate. 

We note that our model calculations used the pre-dimming 
spectrum. It would interesting to see whether the spectrum has 
changed significantly since the dimming. If alumina formed first, 
the peak of the lO^im feature may have shifted to the red. It is 
tempting to associate the alumina with the grey extinction curve, 
but it has little extinction in the opti cal (glass-like ). Better can- 
didates would be Ti02 or pyroxenes ( Woitke 1999). The precise 
agent responsible for the variable polarization also remains unde- 
termined. 



4.3 Evolutionary status 

Jura et al. 12002") pointed out that the luminosity and effective tem- 
perature of L2 Pup do not unambiguously identify it as an AGB 
star, and allow the possibility that it is on the first ascent red gi- 
ant branch. It is not clear whether the spectrum of L2 Pup shows 
technetium, w hich is usually taken as indicating an AGB star. L it- 
tle et al. (1987)^erjorte^a^roteH£_detec^M & 
Hron (1999) listed only a possible detection. 

However, if the if-band has significant obscuration then the 
bolometric magnitude has been underestimated. The bolometric 
correction would also have been overestimated because of the ob- 
served V — K is too large. The two effects combined would mean 
that the actual luminosity is substantially greater than has been as- 
sumed, and is perhaps sufficiently large that placement on the AGB 
is unambiguous. 



Finally, as pointed out by lura et al. (2002 ), the space motion 
of L2 Pup indicates it may be a thick disk star. Given that the thick 
disk ac£comt£ i jOT^H3 i ^o i 235fc - o^ et al. 

2001), it is not unexpected to find a local Mira-like example. If 
L2 Pup is indeed a thick-disk star, as seems plausible but which 
cannot be stated with certainty, the progenitor would likely have 
had a low metallicity and a low mass. 



5 CONCLUSIONS 

Visual photometry of L2 Pup shows an unprecedented dimming 
over the past 5 years. The long-term light curve shows stable pe- 
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Figure 6. Fit to the MIRS spectrum of L2 Pup. The MIRS spectrum (1995) is indicated by the open triangles. Filled squares are the average JHKL 
magn itudes, the IRAS 12 and 25^m fluxes < 1985V the DIRBE 3.5. 4.5. 12 and 25^m fluxes (19901 and the 11.7 and 17.6^m points (2001) of Jura et al. 
(2002). Note that the 12 and 25^tm points refer to very broad bands for which the effective wavelength for this spectral energy distribution is shorter. The left 
plot shows a fit for the silicate feature. The right plot attempts to better fit the longer-wavelength photometry 



riodicity, and we argue that L2 Pup is Mira-like and should be 
classified as SRa. The period stability implies a constant stellar ra- 
dius, which rules out temperature and/or luminosity variations as 
the cause of the dimming. Rather, the dimming seems to arise from 
an episode of dust formation close to the extended atmosphere. 
Episodic dust obscuration events are fairly common in carbon stars 
but have not been seen in (non-symbiotic) oxygen-rich stars. We 
suggest that dust forms continuously but anisotropically, with the 
current dimming event being due to to dust formation along the 
line of sight. 

The red colours indicate reddening from dust, but the extinc- 
tion curve is greyer than found for ISM dust. This could reflect a 
higher fraction of oxides. The change of colour during the dimming 
indicates that already before the dimming, the V-band magnitude 
was significantly affected by circumstellar or atmospheric extinc- 
tion. L2 Pup was one of the few stars located below the Mira P-L 
relation: the derived K-bwad extinction puts the star in closer agree- 
ment with this relation. 

We present a 10-/^m spectrum showing strong silicate emis- 
sion. These observations were carried out in 1995, just prior to the 
recent dimming. The silicate feature can be fitted with a detached, 
thin dust shell, with inner radius 7 x 10 14 cm and outer radius 
~ 1.5 x 10 15 cm. Longer wavelength photometry shows no evi- 
dence for more distant, colder dust. We derive a mass-loss rate of 
M g « 5 x 10 -7 Mq yr _1 , but this value depends on the assumed 
expansion velocity and metallicity — if the dust velocity is high, 
the actual mass-loss rate could be higher. 
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